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The spin- wave excitations of the geometrically frustrated triangular lattice antiferromagnet (TLA) 
CuFe02 have been measured using high resolution inelastic neutron scattering. Antiferromagnetic 
interactions up to third nearest neighbors in the ab plane ( Ji, J2, J3, with J2/J1 ~ 0.44 and J3/J1 ~ 
0.57), as well as out-of-plane coupling (J z , with J z /Ji ~ 0.29) are required to describe the spin- 
wave dispersion relations, indicating a three-dimensional character of the magnetic interactions. Two 
energy dips in the spin-wave dispersion occur at the incommensurate wavevectors associated with 
multiferroic phase, and can be interpreted as dynamic precursors to the magnetoelectric behavior 
in this system. 

PACS numbers: 75.30.Ds, 75.50.Ee, 61.12.-q 



Geometrically frustrated magnetic systems have re- 
ceived considerable attention in recent years due to the 
presence of extraordinary magnetic properties [l], 0, 0, QJ • 
The delafossite CuFe02 is of particular interest because 
of the discovery of multiferroic phenomena with either 
the application of a magnetic field or the substitution 
of Fe 3+ with nonmagnetic Al 3+ ions [1, @, 0, E[. As 
a model material of triangular lattice antiferromagnet 
(TLA), CuFe0 2 forms an Ising-like 4-sublattice (TTil) 
antiferromagnetic order at low temperature, with spin 
moment pointing along the c axis [Figs. 1(b) and 1(c)]. 
This is contrary to other TLA where the three spins align 
at 120° from each other in the basal plane Despite 
the Heisenberg spin of the orbital singlet of Fe 3+ mag- 
netic ions (L = 0, S = 5/2), the metamagnetic tran- 
sition with magnetic field up to 60 T [10[ and the suc- 
cessive thermal-induced magnetic transition are well 
explained with the two-dimensional (2D) Ising model, 
where the Till s P m structure is stabilized via the long 
range magnetic interactions [l2[. Recent results suggest 
that the exchange interaction between adjacent hexagon 
planes in CuFe02 could be equally important as the cou- 
pling within the planes in this seeming 2D system [l3| ]. 
In this Letter, we report high resolution inelastic neutron 
scattering (INS) results of the elementary magnetic exci- 
tations in order to clarify the nature of the magnetic in- 
teractions in this TLA. We found that magnetic exchange 
interactions up to the third nearest neighbor (NN) have 
to be included in order to characterize the dispersion re- 
lations in the basal plane. We also found dispersive exci- 
tations perpendicular to the hexagonal plane indicative 
of a substantial interlayer coupling. Most remarkably, 



we found that the spin-wave dispersion has two dips at 
the wavevectors associated with the multiferroic phase. 
Our observations reveal that the quasi-Ising like spin or- 
der results from the delicate balance between competing 
interactions. The static magnetic order is sensitive to 
external perturbations (magnetic field or impurity), and 
could be transformed into a noncollinear spin structure 
that is intimately related to the polar state. 

A single crystal CuFe02 (mass 4.5 g) was grown by 
the floating zone technique, and additional small crystals 
were crushed to powder for preliminary measurements. 
The INS measurements were performed at the Disk 
Chopper Spectrometer and at the cold neutron triple-axis 
spectrometer SPINS at the NIST Center for Neutron Re- 
search. The Disk Chopper Spectrometer measurements 
on the powder specimen revealed a magnetic excitation 
bandwidth between 0.9 meV and 3.2 meV [Fig. 1(e)]. 
We therefore constrained our triple-axis measurements 
to the energy transfer E < 3.9 meV. Although the crys- 
tal undergoes a lattice distortion from hexagonal R3m to 
a lower monoclinic C2/m symmetry at low T 3, EH, we 
use the hexagonal notation throughout this paper with 
lattice parameter a = b — 3.03 A and c = 17.17 A for the 
sake of simplicity. The crystal was aligned in the (HHL) 
scattering plane. At SPINS vertically focused pyrolytic 
graphite (PG) crystals were used as the monochromator, 
and the final neutron energy was fixed to Ef = 3.7 meV. 
A BeO filter is placed in the scattered beam to remove 
the higher order wavelength contaminations. Horizontal 
collimations of open-80'-sample-80'-open were employed. 

Given the spin structure shown in Figs. 1(b) and 1(c), 
the magnetic Bragg peaks appear at (0.25, 0.25, 3/2) and 
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FIG. 2: (Color online) Representative spin- wave excitation 
scans at T = 1.0 K. Constant wavevector (Q) scans along (a) 
the [if, if, 3/2] direction and (b) the [0,0, L] direction. Solid 
lines are Lorentzian fits. Data are shifted for clarity. 



FIG. 1: (Color online) (a) Schematics of the crystal struc- 
ture and the relevant magnetic exchange interactions between 
neighboring Fe 3+ ions, (b) The ground state spin configura- 
tion in the hexagonal plane. The magnetic unit cell is shown 
by the (yellow) shaded area, (c) The || j j spin structures are 
coupled antiferromagnetically along the c axis, (d) Schematic 
diagram of the (if, K, 0) scattering plane in reciprocal space. 
The structural Bragg peaks are represented in solid symbols, 
the magnetic Bragg peaks displayed in open symbols are pro- 
jections into that scattering plane, (e) Intensity contour map 
of the inelastic measured using the time-of-flight Disk Chop- 
per Spectrometer on a powder specimen (mass ~ 10 g), with 
an incident neutron wavelength of 3.6 A. 



equivalent positions. We first studied the magnetic exci- 
tations within the hexagonal plane along the [if, if , 3/2] 
direction [the solid arrow in Fig. 1(d)]. Figure 2(a) dis- 
plays representative energy scans at constant wavevector 
transfer Qs. The scattering profiles are highly structured. 
In all cases, we identified the lower lying excitation as the 
spin- wave branch in the [if, if, 3/2] direction, while the 
peaks at higher energies originate from unintended scans 
along the open arrow shown in Fig. 1(d) due to the crys- 
tal hexagonal domains [16]. We notice that the dispersion 
does not have a minimum at the magnetic Bragg point 
(0.25,0.25,3/2). Instead, two minima appear at sym- 
metric positions at Qmi = (0.21,0.21,3/2) and Q m 2 = 
(0.29,0.29,3/2), with a energy gap A = 0.9 meV,in 



agreement with the data reported by Terada et al. [17]. 
The scattering profiles along the [0, 0, L] direction, on 
the other hand, show much cleaner features due to the 
absence of twinned domains. The magnetic excitations 
along this direction are clearly dispersive, ranging in en- 
ergy from around 1.3 meV at the zone center to 2.5 meV 
at the zone boundary. This indicates that the previous 
treatment of CuFe02 as a 2D TLA is an oversimplifica- 
tion. 



In order to obtain an expression of the dispersion re- 
lation for CuFe02, we start with the generic Heisenberg 
Hamiltonian: 

H = -\Y. J l & ■ % - \ E 3 to Si-Sj-D^Sl, (1) 



where J2i j indicates summation over pairs of spins, jfj 
is the in-plane exchange interactions (Ji, J 2 , J3 and 
higher order NN coupling as defined in Fig. 1), is the 
interactions between adjacent planes (J z and J' z ) [18[. D 
is the single-ion anisotropy. Using a Holstein-Primakoff 
transformation and a 1/5 expansion, we obtain the ana- 
lytic form of the spin-wave dispersion relation: 

e 2 (k) = A 2 -B 2 ±[(C 2 -C* 2 ) 2 +4\AC-BC*\ 2 } 1 / 2 , (2) 

where 

A = 2S[D - Ji + J 2 (l - cos kyVSa) - J 3 (l + cos2k x a) 

— J z + 2J' Z (1 — cos k z d cos k x a)], 
B = —2S[(J\ + 2 J3 cos k y V3a) cos k x a + J z cos k z d], 
C = -2Scos(k y V3a/2)[J ie iKa/2 + J 2 e- 3iKa/2 } 

-ASJ' z cos k z dcos(k y V3a/2)e- ikxa/2 . 

Here, k x , k y , k z are the wavevectors along the [110], [110], 
and [001] directions respectively, and d = c/3 is the spac- 
ing between adjacent hexagonal planes. 

Along the high symmetry [001] direction, the lower 
branch of the spin-wave dispersion relation described in 
Eqn. (2) reduces to the much simpler form: 



e 2 (k z ) 



4S 2 [(D - G)(D - G + F) - 
FG cos k z d — G 2 cos 2 k z d] , 



(3) 



where F = 2(J 2 - 2J 3 ) and G = J z - 2J' Z . The band- 
width of e 2 (k z ) is mainly controlled by the out-of-plane 



3 



(q.q.o) (riu) 




H (rlu) 



FIG. 3: (Color online) (a) Spin-wave dispersion along the 
[H,H, 3/2] direction. Points are experimental data and the 
solid line denotes a fit to the model described in the text. Only 
the data from the lowest branch are included in the analysis. 
The error bars represent one standard deviation, (b) Pre- 
dicted dispersion relations along the open arrow in Fig. 1(d) 
(contribution from the twinning domains), (c) Phase dia- 
gram from Ref. [l2|- The solid symbol is the value obtained 
from our measurements, (d) The structural super lattice peaks 
appear at Ql = 2Q rn and Ql — 1/2 — 2Q rn (due to ex- 
change striction) when the incommensurate spin order occurs 
at /iqH — 9 Tesla. 



interactions J z and J' z [i.e. £ 2 (1.5) -£ 2 (3.0) = -85 2 FG], 
while the single-ion anisotropy term D is determined by 
fitting the experimental data to the reduced dispersion 
relation Eq. (3). We obtained DS = 0.17 meV, which 
is larger than the estimated value of DS = 0.05 meV 
from magnetization measurement |l3[. The solid lines 
in Fig. 3(a) and Fig. 4 represent the global fits to the 
observed spin- wave energies along [H,H, 3/2], [0,0, L] 
and [q,q,L] (with q = 0.21 and 0.25) using the model 
described above. Long-range antiferromagnetic interac- 
tions (up to Js) are found (Ji5 = —1.14 meV, J2S = 
—0.50 meV, and J3S = —0.65 meV) to satisfactorily de- 
scribe the in-plane dispersion relation [19[. The inclu- 
sion of higher order interactions (J4S) in our analysis 
did not improve the fits. Only one NN exchange inter- 
action (J Z S = —0.33 meV) is necessary to characterize 
the spin- wave dispersion out of the basal plane, the ad- 
dition of a second out-of-plane exchange interaction also 
yields a small value (J' Z S = 0.01 meV) with no significant 
improvement of the fit. 

We now compare our fitting parameters with the ear- 
lier theoretical calculations. The phase diagram of the 




FIG. 4: (Color online) Spin-wave dispersions along the (a) 
[0,0, L], (b) [0.21, 0.21, L], and (c) [0.25, 0.25, L] direction. 
Points are experimental data and the solid line are fits to 
the model described in the text. Only the data at the lowest 
branch are presented in (b) and (c). 



2D Ising spin TLA has been investigated by Takagi and 
Mekata using Monte Carlo simulations [12]. They found 
that the collinear 4-sublattice structure is only stable 
with fairly large J 2 and J3 contributions. The phase 
is bound within a narrow range of —0.5 < J2/IJ1I < 
and J3/J2 < 0.5 [the shaded area in Fig. 3(c)]. Our 
result, with J2/J1 =0.44 and Js/ Ji=0.57, is positioned 
within this shaded region near the boundary where the 4- 
sublattice state coexists with the 8-sublattice state [l2j. 
The phase diagram is not affected by the inclusion of the 
single-ion anisotropy D, as this term uniformly shifts the 
energy at each spin site. Similarly, this phase diagram 
remains unchanged with respect to the NN out-of-plane 
exchange coupling J z , since the total energy stays invari- 
ant regardless of the detailed in-plane magnetic structure, 
providing there is a negligible J' z . The agreement be- 
tween our estimation of the in-plane exchange constants 
and those predicted by Takagi and Mekata implies that 
the mean field approach can effectively explain the emer- 
gence of the long-range magnetic order. 

Another important result is the observation of two 
energy dips at Qmi an d Qm2 in symmetric positions 
around the commensurate magnetic Bragg point r = 
(0.25,0.25,3/2) characteristic of the 4-sublattice phase. 
Remarkably these are precisely the same wavevectors 
associated with the incommensurate magnetic order of 
CuFe02 under a magnetic field [20[, and are dynamic 
precursors to the multiferroic phase. By decreasing the 
single-ion anisotropy D without modifying the exchange 
constant JVs, one can essentially close the spin gap A, 
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while the characteristic Q m related to the energy min- 
imum remain unchanged. In any antiferromagnet, the 
spin- wave modes are linearly split in a magnetic field and 
a critical value of H C (T) = A/gfiB will destroy the local 
stability of the antiferromagnetic state. For A = 0.9 meV 
this critical value is 7.7 Tesla, which is in agreement 
with the experimental observation that incommensurate 
magnetic order appears for fi^H > 7.0 Tesla [5]. The 
softening of the magnetic dynamics achieved by apply- 
ing a magnetic field reveals a close connection between 
the breakdown of the Till spin structure and the forma- 
tion of a new type of magnetic order. The local stability 
of the antiferromagnetic state can also be destroyed by 
chemical substitution. A polar state has been recently 
discovered in the doped CuFei_ x Al x 02 in zero magnetic 
field possessing an incommensurate spin structure with 
the same Q rn [3, Q . 

There remain open questions in this TLA system. For 
instance, why are the exchange interactions up to J3 re- 
quired to stabilize the collinear spin structure? Although 
the slight increase in bonding angle of Fe-O-Fe along 
the [110] direction [Fig. 1(b)] below the Neel tempera- 
ture might explain a dominant antiferromagnetic inter- 
action Ji, there is no obvious justification for the exis- 
tence of comparable AF interactions between higher or- 
der neighboring spins. Second, the microscopic origin of 
the anisotropy gap A is yet to be understood. One would 
expect the magnetic properties should be isotropic since 
Fe 3+ is an orbital singlet. One possible explanation might 
be related to the lattice distortion accompanied with the 
magnetic order which induces an anisotropy between the 
Fe-3d and 0-2p hybridization as well as the charge trans- 
fer between them [21]. Finally, the excitation data are 
not completely characterized by the Heisenberg model 
described in Eq. (1). Figure 3(b) shows the expected dis- 
persion curves from the twinned magnetic domains. Only 
some features agree with the experimental data. Asym- 
metric dispersion curves appear with local minima lo- 
cated at Q = (0.22,0.22,3/2) and (0.36,0.36,3/2). How- 
ever, the magnetic gap and overall energy bandwidth do 
not agree with the observations. 

In summary, the spin- wave spectrum of the geometri- 
cally frustrated TLA CuFe02 has been mapped out us- 
ing INS. We observed a considerably dispersive excitation 
along the c axis, which reflects the 3D character of the 
magnetic interactions. We have also determined the rel- 
evant spin Hamiltonian parameters, confirming that the 
next NN and the third NN interactions are required to 
explain the spin dynamics. Finally, the local minimum 
of the dispersion curves reveal the dynamic precursory to 
the multiferroic phase. 
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